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We present a theory of the quasi two-dimensional electron gas (2DEG) systems which appear near 
the surface of SrTiOa when a large external electric field attracts carriers to the surface. We find 
that non-linear and non-local screening by the strongly polarizable SrTiC>3 lattice plays an essential 
role in determining 2DEG properties. The electronic structure always includes weakly bound bulk- 
like bands that extend over many SrTiOa layers. At 2D carrier-densities exceeding ~ 10 14 cm -2 
tightly bound bands emerge that are confined within a few layers of the surface. 

PACS numbers: 73.20.-r, 71.10.Ca, 68.47.Gh, 74.78.-w 



I. INTRODUCTION 

Two-dimensional electron gases can be formed in 
SrTiC>3 crystalline by gatin gP^ b y forming an inter- 
face with a polar perovskitej 1 * 3 * 7 * 8 * or by placing a 6- 
doped 9 layer inside a bulk crystal. Although 2D elec- 
tronic systems at LaAlOs/SrTiOs interfaces have re- 
ceived particular attentionpED there has also been im- 
portant progress with other material systems™ SrTiC>3 
two-dimensional electron gases (2DEGs) appear to be 
strongly correlated when their thermodynamics is probed 
capacitiveljW and exhibit both superconductivity^ and 
magnetism, 1 ^ sometimes simultaneously.^ There is at 
present only a very primitive understanding of the mea- 
sured properties of these potentially interesting 2DEG 
systems. The current paper is motivated by the view 
that progress can be accelerated by the development of 
concrete microscopic models that are simplified relat ive 
to full ab-initio electronic structure calculations J1SUSI al- 
lowing electric properties to be estimated easily and com- 
pared with experiment. 

In this paper we present a model of SrTiC>3 2DEGs 
that is partly phenomenological and simplified, but still 
sufficiently realistic to be predictive. We focus on electro- 
statically gated surface 2DEGs, although our approach 
applies without much change to the case of interface- 
confined systems. The same model is readily adapted 
to describe <5-doped 2DEGs inside the STO bulk, STO 
2DEGs that are modulated by a back gate, and 2DEGs 
in other d? systems, for example KTa03. The model 
assumes that the itinerant electronic degrees of freedom 
are derived from the SrTi03 t2 g bands. We use a near- 
est neighbor tight-binding model to describe hopping 
between Ti02 planes and either tight-binding or k ■ p 
models to describe wavefunction variation within Ti02 
planes. The strength of inter-plane hopping parame- 
ters, and the values of the heavy and light masses within 
planes are estimated on the basis of recent ARPES^ and 
bulk magnetic oscillation^ experimental results. Some 
aspects of the 2DEG electronic structure are sensitive to 
the influences of spin-orbit coupling and SrTiOs's low- 
temperature tetragonal distortion on the host material's 
conduction band, even at the highest 2D carrier densities. 



The extremely strong dielectric response of the SrTi03 
lattice plays a key role in our model at all carrier densi- 
ties. Our main results are summarized in Fig.[T] We con- 
clude that unless vertically confined on both sides by vac- 
uum or insulating tunnel barriers, SrTi03 2DEGs spread 
across a large number of Ti02 planes. This property is a 
direct consequence of the host material's very large lin- 
ear dielectric constant, which weakens confinement, and 
occurs in spite of relatively large carrier masses which 
favor confinement. At high carrier densities, and hence 
large electric fields, dielectric screening saturates and the 
2DEG is mostly confined to the first few Ti02 planes. 
However a portion of the 2DEG, making a contribution 
to the 2D density that is approximately fixed in abso- 
lute terms, still spills over many layers. This low-density 
weakly confined part of the 2DEG can make an impor- 
tant contribution to some 2DEG properties. 

Our paper is organized as follows. In the following 
section we provide a detailed explanation of the model 
that we use. We have identified three different density 
regimes for SrTi03 surface-bound 2DEGs. In Sees. III-V 
we characterize the nature of the 2DEG electronic struc- 
ture in low (n < 1 X 10 14 cm~ 2 ), mid (1 x 10 14 cm~ 2 < 
n<5x 10 14 cm- 2 ), and high (n > 5 x 10 14 c m - 2 ) 2DEG 
carrier density regimes respectively. Finally in Sec. VI 
we summarize our results and speculate on the types of 
electronic properties that might be achievable in SrTi03 
2DEG systems. 



II. MODEL 
A. t2 g tight-binding model 

SrTi03 is a non-polar pseudo-cubic band insulator 
with an electronic gap of ~ 3.2 eV separating its oxygen 
p-orbital dominated valance band, from its Ti i2g-orbital 
dominated conduction band. The Ti e g bands are split 
by the crystal field and pushed up in energy by ~ 2 eW^* 
relative to the ti g bands and are therefore neglected in 
our model. The conduction band minimum is at the T 
point. The bulk t2 g bands are split at the T point, in 
the first place by spin-orbit interactions which push one 
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FIG. 1: (Color online) Orbital character at k± — of occu- 
pied doubly degenerate 2D subbands for a series of total areal 
densities. The xy and {yz,zx} fractions in the orbital char- 
acter of each band are represented by light blue and light red 
shading respectively and the band indices are ordered from 
lowest energy to highest. Both spin-orbit interactions and 
tetragonal splitting have been included in the band-structure 
model. The percentage of the total density associated with a 
given subband is summarized in the inset piechart, in which 
the rings from inside to outside correspond to the density val- 
ues from lowest to highest. Individual band contributions for 
a particular density are ordered lowest energy to highest in 
a clockwise direction. The red lines in the pie chart separate 
the low, mid, and high density regimes identified in the text. 
At the highest total densities, most electrons occupy a small 
number of strongly confined bands. 



time-reversed pair of ti g bands up by ~ 18 meV relative 



to the other two pairs. The two lower energy bands are 
further split below 110 K by a rotation of the octahedral 
oxygen environment which is responsible for a tetragonal 
distortion.^ 

The starting point for any phenomenological descrip- 
tion of the electronic structure of SrTiC>3 2DEGs is an 
accurate representation of the bulk material electronic 
structure. Although this material has been studied for 
many years, its conduction band minimum was until re- 
cently not characterized with an accuracy sufficient to 
model low-carrier density 2DEGs. To remedy this, Allen 
et al. conducted^ magnetotransport studies on a se- 
ries of low density electron-doped MBE-grown samples 
of SrTiC>3 in fields up to 31T and fit a 6 band k ■ p model 
of the Ti t2 g bands to the magnetic oscillation data. The 
bulk band parameters used here have been taken from 
that study. These experiments suggest that in the ab- 
sence of spin-orbit coupling the tetragonal distortion at 
4 K would push the xy bands up by ~ 3.2 meV rela- 
tive to yz and zx bands. Together the two-corrections 
fully lift the ti g manifold degeneracy, even in the bulk. 
(Because of orbital mixing, the spin-orbit (SO) splitting 
and tetragonal splitting parameters do not directly cor- 
respond to the r point band energies.) 

Although the p—d oxygen bonding and e g anti-bonding 
orbitals do not explicitly enter our model, they do appear 
implicitly in the form of the Hamiltonian. Consider for 
example hopping between a Ti xy orbital and its neigh- 
boring oxygens within a Ti02 plane (Fig. [2^,). Along 
the x-direction, the dominant bonding is tt pc i through the 
O-py orbital and along the y-direction, 7r pt ; through the 
O-pa; orbital. Overlap with other O-p orbitals is small by 
symmetry. This rule is preserved througout the Ti - O - 
Ti bonding network. For the Ti - yz orbital, ir p d bonding 
dominates along the y-direction through 0-p 2 orbitals. 
Bonding along the x-direction vanishes in a Slater-Koster 
two-center approximation and is weak. Temporarily ig- 
noring the spin-orbit and tetragonal distortion effects, 
these observations suggests a tight binding model for a 
single isolated layer of the form 



-2t'cos{k x ) — 2tcos(k y ) \ ( yz, a 

H S G L = ( -2tcos(k x )-2t'cos(ky) W zx,a )■ , (1) 

—2tcos(k x ) — 2tcos(k y ) I I xy, a 



where the cubic lattice constant is used as a length unit, the metal lattice site energies are used as the zero-of-energy, 
t quantifies the dominant ir p d bonding process and t' describes the weaker bonding process. The column on the right 
specifies the orbital representation used for this Hamiltonian matrix. Hopping terms that couple different <2g orbitals 
are allowed^ from a symmetry point of view. However, Allen et al. were unable to distinguish this mixing parameter 
from zero in their recent analysis of SdH data!^ We therefore ignore these processes in our model. For lower carrier 
densities it is sometimes convenient to use a simplified version of this model in which we expand Eq. [I] around the 
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2D r point. We find that for 2D wavevectors that are small compared to Brillouin-zone dimensions 

( €yzfl+t'kl + tk 2 y \ r yz,<J ) 

H? L = I e ZXt0 +tk 2 x +t'k 2 y zi,(7 , (2) 

V e xyfi + tk 2 x + thy J I xy, a I 




where e VZi o — e zx ,o = —2t — 2t' and t X y,o = — 4i. We 
use this low density form for the planar Hamiltonian for 
most of the calculations presented below. The more gen- 
eral tight-binding model must be used when 2D carrier 
densities are large and confinement is strong, and can be 
adopted when required without essential complication. 

In the same representation, adjacent 2D layers are cou- 
pled by an interlayer hopping term of the form 



(3) 



Here the symmetry of the bonding network has again 
been employed to note that the xy orbital has the weaker 
interlayer coupling, t' . Because t' is expected to be sub- 
stantially smaller than t, the xy bands in single-layer 
2DEGs are pulled down by ~ 2t at the T point relative 
to the yz and zx bands. In the bulk limit, on the other 
hand, the three bands are degenerate because each has 
two strong hopping and one weak hopping direction. Any 
amount of confinement in the z direction pushes the bot- 
tom of the {yz, zx} bands up relative to the xy band and 
leads to orbital polarization. 

For low-carrier densities on-site (fc-independent) terms 
due to tetragonal distortions and spin-orbit coupling 
must be included.^ The tetragaonal distortion is repre- 
sented by a parameter Ap which characterizes the differ- 
ence in site energy between xy and {yz, zxj-orbitals, and 
spin-orbit coupling by an interaction strength parameter 
A50 The distortion Hamiltonian is 





A T 




(4) 



and the spin-orbit Hamiltonian, modeled in an atomic 
approximation, is 
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(5) 



Our model for the electronic structure of SrTiO,3 
2DEGs combines the single-particle model explained 
above with a Hartree approximation for electron-electron 




-Qui -Q(u*2-u1) 

FIG. 2: (Color online) (a) p — d bonding network for Ti- 
ll/ orbital. The dominant bonding is in-plane 7r p d while the 
out-of-plane bonding is weak, (b) Schematic representation 
of the electric field drop in our model due to lattice screening 
(dashed blue) and electronic screening (solid red). Numerical 
factors have been dropped. 



interactions. The external electric field which produces 
surface confinement is screened by carriers and by lat- 
tice relaxations of the partially ionic SrTi03 crystal. In 
SrTi03 lattice screening is strong and non-linear and 
plays a subtle and essential role in confinement. 



B. Lattice relaxation model 

The exceptionally strong and temperature-dependent 
linear dielectric response of bulk SrTiOs is associated 
with a soft optical phonon mode in which positively 
charged Sr and Ti atoms move in opposition to the 
negatively charged oxygen octahedra. Displacement of 
this mode in response to an external electric field pro- 
duces screening. Because the mode is extremely soft 
only near the center of the Brillouin-zoneJ22l it responds 
strongly only when a large external field persists over sev- 
eral Ti02 layers. In addition this screening response is 
very non-linear, saturating at very large electric fields. 
Since the reduction in electric field is proportional to 
the phonon mode displacement, saturation occurs be- 
cause the phonon mode is anomalously soft only for small 
displacements f22l In an attempt to capture this behavior 
qualitatively, we use a simple model of lattice relaxation 
which focuses on the soft-mode only. We write the lattice 
energy as 



U 



(6) 
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where itj is the displacement of the soft-mode on the 
i-th lattice site, -E 1 , is the average electric field in the 
i th cell, Q is an effective charge defined in terms of 
the polarization-density per unit soft-mode displacement, 
and 7 is a parameters chosen to capture the non-linearity 
of the dielectric response as discussed further below. Here 
Ki j = Ki_j is the dynamical matrix at 2D wavevector 
q± = 0. We fit to the soft-mode phonon disper- 
sion using a form with a local on-site contribution and 
a Gaussian non-local contribution. In momentum space, 
this takes the form, 



r2 _ f 2-^(q+G) 2 _ f 2-^(q+G) 2 



JO 



Jl< 



J2> 



K(q,G) = (27r) 2 M| 

(7) 

where q is the lattice momentum, G is a reciprocal lattice 
vector, and fi = 24amu is the appropriate reduced mass 
for the Ti atom moving opposite the oxygen octahcdra. 
The parameters /o , the strength of the on-site term, and 
fx and ot\ were chosen to reproduce the measured phonon 
dispersion. Because no low temperature phonon data 
exists in the literature, f% and a.2 have been added to 
capture the low temperature dielectric response of the 
bulk material. By minimizing Eq. [6] in the absence of an 
electric field and evaluating Eq. [7] at q = we find that 
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(8) 



where u^l is the mode displacement at which non- linear 
dielectric response is seen. (See the discussion below). 

Because the relative displacements of all atoms are 
known, only a single displacement coordinate, it;, is 
needed to describe the response of the unit cell to per- 
turbations along the principle crystal axes. Given the 
electric field in each cell of the crystal, Eq. [6] can be min- 
imized to find the appropriate set of displacements, {it*}- 
We define the three-dimensional polarization density of 
the SrTi0 3 as 



Pi = ~^Qu* [e(zi - a/2) - e( Zl + a/2)} (9) 
a* 

where Zi is the location of the Ti02 layer of interest, a is 
the lattice constant of the crystal, and &(z) in the Hcav- 
iside function. The precise way in which the polarization 
density is mapped onto our lattice model is immaterial 
on length scales larger than a lattice constant. To find 
the effective charge parameter we use the standard def- 
inition of the screened electric field and linear dielectric 
constant, 



E = Eq + AttP « eE . 



(10) 



After minimizing Eq. [6]in the linear, bulk limit and using 
the definition of the polarization from Eq. [9j we find that, 



4-7T 



(6-1). 



(11) 



To make contact with the measured properties of the bulk 
material in a straightforward way, we use 90K values 
for the phonon dispersion^ and dielectric constant.^ In 
terms of model parameters ljx — 2tt(/o — fx) and the 90K 
dielectric constant are given in Table 1. With this we find 
Q = 8.33e, a value comparable to those used in models 
of this type for bulk SrTi0 3 P 

The electric field in Eq. [6] can be found by solving 



V -E{z) 



Aire 



^rii5(z-Zi)+An^V -Pi, (12) 



with the boundary conditions that E{— oo) = Eq and 
E^ = 0. The electric field boundary conditions are dis- 
cussed below. In Eq. 12 e is the electron charge, f M is 
the high frequency dielectric constant due to electronic 
screening, and rij is the number density of itinerant elec- 
trons in Ti02 layer i. Both lattice relaxation and con- 
duction band charge accumulation screen the external 
electric field. This is represented pictorially in Fig. [2)3. 



C. Electric Field Boundary Conditions 

In the calculations presented below we assume that the 
electric field above the surface of the SrTi0 3 has been 
set experimentally either by gating or by forming an in- 
terface with a polar dielectric.^ In the latter case Eq 
is ideally set by the polarity of the material, but can 
also be influenced by surface reconstructions or other de- 
tailed material considerations that can be sensitive to 
uncontrolled aspects of growth. Because we have gated 
systems in mind, we consider that Eq can be varied ex- 
perimentally over a broad range. In this calculation we 
set the electric field below the SrTi0 3 2DEG, Eb u ik, to 
zero, assuming that the sample lies on a grounded metal- 
lic substrate. (If the SrTi0 3 sample was thin, a con- 
ducting layer under the sample could be used as a gate 
and Ebuik could be varied.) By integrating the Poisson 



equation (Eq. 12 ) and noting that the lattice relaxation 
contribution to E{z) vanishes far below the surface when 
Ebuik —> 0, we conclude that the 2DEG density is fixed 
by Eq alone: tit — n i — CooEq/ Aire. We can therefore 
replace this parameter by the total 2DEG density tit and 
present results as a function of that parameter. 

We incorporate the layer-dependent electric poten- 
tial contribution to the Hamiltonian by integrating E(z) 
across the 2DEG to obtain a layer-dependent potential 
Vi which must be determined self-consistently along with 
the 2DEG density-distribution and the soft-mode dis- 
placement field. With this, the Hamiltonian of the sys- 
tem becomes, 
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Model Parameters 


Lattice Constant 


a 


3.904 A 




t 


236meV 


Electronic Parameters 


t' 

A so 


35meV 
18meV 




A T 


3.2meV 




eo 


24408 




£l 


1340 


Dielectric Response 


Q 


5.5 
8.33e 




UNL 


0.0034A 




fo 


4 x 10 i2 c/s 




fl 


2.73 x 10 12 c/s 


Dynamical Matrix 


h 


0.97 x 10 12 c/s 




«1 


1.15a 




«2 


5a 



TABLE I: Parameters used in the current study. The elec- 
tronic structure parameters have been taken from Ref. 1181 
while eoo, ei and eo were taken from Refs. I22l fc l24l 



H = 



<i,J> 



H 



SL 



T + H 



so 



Vi). 



where the double sum in the first term is over neighboring 
layers. In Eq. |13[ we work in the representation c = 
{ c xy,f> c yz,^7 c zx _-f, c-xy,\.i c yz,i, c zx,i} s ° that H so has the 
form of Eq. [51 

The layer resolved density rij =< c-c,; > is calculated 
from Eqs. ( |12|13 1 and minimization of Eq. [6j We have 
carried these self-consistent field calculations to conver- 
gence over a wide range of densities for a system that is 60 
unit cells wide. Because of the long tail in the density- 
distribution discussed at length below, it is difficult to 
achieve self-consistency and we were forced to mix in no 
more than ~ 1% of new results in the iterative update 
of the charge density. Although the model described in 
this section is crude in some respects, certainly crude 
compared to ab initio electronic structure calculations 
with full lattice relaxation, it is strongly motivated by 
the cumbersome character of the fully microscopic calcu- 
lations under these circumstances. The model could be 
made more quantitative by being bench-marked against 
ab initio calculations or, perhaps more reliably, by com- 
parison with experiment. 



III. LOW CARRIER DENSITIES 

1 x 10 14 cm -2 < tit 



For the circumstance considered here the total carrier 
density is proportional to the electric field just above the 
SrTi03 surface and the largest internal electric fields are 
closest to the surface. We define the low-density regime 
by the requirement that the largest electric fields are 
smaller than the scale at which non-linear screening be- 
comes important. This field scale is set by the model 




FIG. 3: (Color online) (a) Fermi surface and 2D band struc- 
ture for 72t = 8.3 x 10 12 . The dominant orbital character at 
the 2D F point is represented for each band by the color of 
the bands and Fermi lines - blue and red for xy and {yz, zx}, 
respectively. The zero of energy is set to the potential min- 
imum in the first layer. The Fermi energy is indicated by a 
solid (black) horizontal line. 



parameters ujvl and eoo, which can be determined ap- 
proximately from experiment. We have estimated ujvl 
by comparing Eq. [9] with the deviation from linear re- 
sponse seen in the polarization of bulk STO crystals.^ 
This value is listed along with other model parameters in 
Table I. The model parameters we have chosen reflect the 
estimate that non-linear screening becomes important for 
carrier densities larger than ~lx 10 14 cto~ 2 . 

For linear screening some 2DEG properties are similar 
to those of covalent semiconductor 2DEGs and can be 
estimated following the same lines as in Stern's pioneer- 
ing study of Si/Si0 2 MOSFET 2DEGsP In particular 
the confinement length scale w can be crudely estimated 
by equating the quantum confinement kinetic energy cost 
and the confinement electric potential scale. Neglecting 
numerical factors we therefore set 



cEqw Aire 2 nTW 



(14) 
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to obtain 
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4Ttme 2 riT 



1/3 



(15) 



In the linear screening regime the confinement length 
scale decreases quite slowly with the total 2DEG den- 
sity. The hopping parameters of Table I can be con- 
verted to effective masses for the t 2g bands; the light 
mass that describes the vertical confinement of the most 
poorly confined {yz, zx} bands is ~ mo where mo is the 
bare electron mass. When combined with the extremely 
large low-temperature bulk dielectric constant of SrTiC>3 
(e ~ 25000), we estimate that w is close ~ 50 SrTiC>3 
unit cells even at the top end of the low-density regime. 
We therefore expect that the hard wall at 60 unit cells 
used in our calculations influences our numerical results. 
The main point of these qualitative considerations is that 
we should expect weak surface confinement at low carrier 
densities because of very strong dielectric screening. 

In Figs. [| we illustrate a typical 2D band structure in 
the low density regime. Here the bottom band is begin- 
ning to reflect the increase in xy character expected from 
confinement, and the small size of the subband splittings 
is in qualitative agreement with the estimated scale of 
size-quantization energies: 



• Aso = 0, At = 
■ Aso * 0, At = 

♦ Aso * 0, At * 




lxlO 12 



lxlO 13 lxlO 14 
Density [cm -2 ] 



high 



10~ 4 eV. 



mur 



(16) 



The small subband splittings imply that the 2DEG is 
3D in character unless temperatures are low and disor- 
der extremely weak. The vertical spread of the 2DEG is 
expected to get smaller, and the subband splitting larger 
with increasing temperature as the dielectric constant 
value decreases 

Low carrier density properties are strongly influenced 
by spin-orbit coupling which hybridizes the t2 g basis 
states and induces a splitting at the T point in the bulk. 
One effect of spin-orbit coupling is to weaken the 2DEG 
surface confinement by hybridizing xy bands with yz and 
zx bands that have smaller masses perpendicular to the 
surface. Spin-orbit induced hybridization allows the xy 
subbands to communicate between layers through their 
{yz, zx} admixtures which are less easily confined. In 
the low density regime, spin-orbit splitting is pronounced 
enough to change the dominant orbital character of the 
2D subbands. 

The tetragonal distortion increases the site energy of 
the xy band - further enhancing the role of the less 
confined {yz, zx} components. (We have assumed that 
the tetragonal axis is parallel to the surface normal, as 
expected near a surface.) Initially, confinement energy 
scales are weak compared to the tetragonal splitting en- 
ergy. As the carrier density and the energetic width of 
the occupied states increase, spin-orbit coupling becomes 
less important and the xy fraction of the lowest energy 
most highly occupied band increases. (See Fig. [I]) The 
influence of the spin-orbit and tetragonal splittings on the 
spatial distribution of electrons is summarized in Fig. [4] 



FIG. 4: (Color online) Average separation from surface, (z) 
and standard deviation a of the electron distribution across 
layers as a function of total density, and its dependence on SO 
splitting and the tetragonal distortion. The tendency of SO 
coupling and tetragonal splitting to weaken surface confine- 
ment is suppressed when densities reach the mid range. When 
confinement energy scales are not strong enough to overcome 
the tetragonal distortion, SO and tetragonal-splitting induced 
hybridization decrease the spread of the {yz,zx} bands (see 
text). 




Ixl0 l2 2xl0 12 5xl0 12 lxl0"2xl0 13 5xl0 13 

n2D[cnr 2 ] 

FIG. 5: (Color online) Calculated 3D density in the low den- 
sity regime. The relationship between 2D density and average 
3D density follows a 4/3 power law that is consistent with Eq. 
1141 is correct. 
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FIG. 6: (Color online) Lattice displacement as a function 
of total 2D density. For each density, the lattice displace- 
ment of layer i (counting from the surface) is plotted as the 
height of the i th bar segment (counting from the bottom). 
For low-densities, the lattice displacements are small and in 
the linear regime. For mid-range densities, lattice displace- 
ments are suppressed by non-linear screening effects near the 
surface. Weaker lattice screening results in stronger confine- 
ment, larger 2D subband separations, and fewer occupied 2D 
subbands. 



Estimating 3D densities using n^u = ri2Dl{z) where (z) 
is taken from Fig. [4] and n2D is the total density in the 
linear screening spatial region, we find a power law of 
4/3. This suggests that the qualitative estimate of Eq. 
|14|is accurate when screening is linear. 



IV. MID-RANGE DENSITIES: 

1 x 10 14 cm -2 < n < 5 x 10 14 cm 2 

We define the mid-range of densities as that in which 
lattice screening is markedly reduced because of non- 
linear dielectric screening (see Fig. [6j) Because the elec- 
tric field is larger closer to the surface, non-linear screen- 
ing is more important there. The strong surface electric 
fields cause a large fraction of the total electron density 
be confined close to the surface, and size-quantization 
effects to increase much more rapidly with carrier den- 
sity than would be suggested by Eq. [l4j Even though a 
substantial fraction of the total charge density starts to 
become confined within the top few layers, there is still a 
wide tail in the density distribution in the spatial region 
over which the external electric field has been reduced 
to a value less than ~ 10 14 cm _2 e oo /(47re) so that the 
screening is locally linear. In our model the non-locality 
of these screening properties is set by the width in mo- 
mentum space of the long-wavevector limit of the soft 
mode. In our numerical calculations, this low-density 
quasi 3D regime is influenced by our hard- wall cut-off of 



^ 



^180 
— 140 
^ 1 00 

CD 

lS 60 




k x [l/a] 



FIG. 7: (Color online) (a) Fermi surface and 2D band struc- 
ture for tit = 2.0 x 10 14 . The dominant orbital character of a 
band at the 2D F point is represented by line color with blue 
and red indicating xy and {yz,zx}, respectively. The zero of 
energy is set to the potential minimum in the first layer. The 
Fermi energy is represented by a solid (black) horizontal line. 
Although the separation between the lowest energy 2D sub- 
bands is large, many low density subbands with small energy 
separations are still present near the Fermi energy. 



the 2DEG beyond a width of 60 unit cells. 

As was the case in the low density regime, the inclu- 
sion of SO coupling and the tetragonal distortion alters 
the the orbital character of the lowest energy band and 
decreases its surface confinement. Their influence is re- 
duced compared to the low-density regime however. As 
illustrated in Fig. [7] we find that, at the 2D T point, 
the two lowest bands are dominantly xy in character and 
that the next occupied subbands are {yz, zx} in charac- 
ter. Although the number of 2D subbands has increased 
significantly, only a few are needed to account for the 
most strongly confined part of the density (see the inset 
of Fig. [TJ. 
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FIG. 8: (Color online) Layer resolved density (as a percent of 
total density) as a function of the total density. For each den- 
sity each segment represents the percent of the total density 
in the corresponding Ti02 layer, starting with the first layer 
and moving upward. For low densities the 2DEG is spread 
over many layers. Above tit ~ 10 14 cm -2 , the confinement 
becomes pronounced. In the high density regime, more than 
50% of the total density is confined within the first layer due 
mainly to reduced lattice screening at large electric fields. 
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The 2DEG electronic structure simplifies again in the 
high-density limit, which we define as the limit with more 
than half of the total density in the first TiC>2 layer - see 
Fig. [8] For large electric fields, and therefore large carrier 
densities, lattice screening is irrelevant near the surface. 
Because of the relatively large conduction band masses, 
compared to typical covalent semiconductor cases, and 
the much stronger electric fields at these carrier densities, 
surface confinement occurs on an atomic length scale. 
The T-point splitting becomes comparable to the single 
layer limit of Eq. [2] - see Fig. [9j While the SO coupling 
leads to hybridization and a decrease in confinement that 
mainly affects the quasi 3D tail of the electronic distri- 
bution, the tetragonal distortion does not have a notice- 
able effect. In this regime, if we neglect the quasi 3D 
tail region, there are only a few spin-degenerate 2D sub- 
bands contributing to the density. At the 2D V point, the 
first three bands are dominantly xy and {yz, zx} - going 
from low to high energy. For the high-density regime, the 
tight-binding model of Eq. [I] must be used. 



VI. SUMMARY AND DISCUSSION 

Using a simplified tight-binding model for the t 2g 
bands, we find that non-linear screening plays an es- 
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FIG. 9: (Color online) (a) Fermi surface and 2D band struc- 
ture for tit = 5.9 x I0 14 . The dominant orbital character of 
the band at the 2D T point is represented by line color with 
blue and red for xy and {yz,zx}, respectively. The zero of 
energy is set to the potential minimum in the first layer. The 
Fermi energy is represented by a solid (black) horizontal line. 
The 2D subband splitting is ~ 200meV - becoming compara- 
ble to the single layer model of Eq. [JJ Near the Fermi energy, 
many subbands are present with ~ meV splitting, which con- 
tribute to the carrier density far from the surface. 



sential role in determining the electronic properties of 
surface confined 2DEGs in grounded SrTiOa. For low- 
density (tit < 10 14 cm~ 2 ) 2DEGs, electrons are dis- 
tributed over many layers because surface confinement 
is weakened by the host material's extremely large linear 
dielectric constant. In the mid and high density regimes, 
a low density tail is still present over 50 or more layers but 
a high density region emerges over the first few atomic 
layers. Although many 2D subbands are always present 
in the mid and high density regimes, the few lowest en- 
ergy subbands which contain strongly confined orbitals 
account for most of the total density (see Fig. [lj. Sub- 
bands that have substantial {yz, zx} orbital character are 
present at all densities in the grounded configuration in- 
vestigated in this study. The presence of many subbands 
with different orbital character at the Fermi level sug- 
gests that the interpretation of transport properties in 
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these systems is not likely to be straightforward. 

SrTiOa is well known for exhibiting superconductivity 
in the bullP^ where it appears over a broad range of car- 
rier densities from ~ 10 18 cm -2 to more than ~ 10 20 cm -2 
and has a maximum value ~ 400mK!2£' In one study of 
gated SrTiC>3 2DEG systems,^ with Hall densities in the 
range from ~ 10 13 cm~ 2 to ~ 10 14 cm~ 2 , the supercon- 
ducting transition temperature initially increased with 
carrier density but decreased at the higher densities. Ref- 
erencing to Fig. [5] for 3D densities associated with the 
weakly confined tail, we find that the measured surface 
2DEG T c 's compare well with values reported for bulk 
systems.^ In another studyP superconductivity was seen 
only at a Hall density of 3.9 x 10 14 cm~ 2 with a transition 
temperature of near ~ AOOmK. Because of experimen- 
tal limitations superconductivity was not seen at other 
densities but could have been present at lower transi- 
tion temperatures!^ In LaAlOs/SrTiOa systems the re- 
ported 2DEG densities are in the low density range. It 
is therefore not surprising that the measured T c values 
are correspondingly suppressed! 12 * 14 ^ (The carriers found 
in LaA103/SrTi03 2DEG systems are thought to be in- 
duced by a polarization discontinuity, although the small 
value of the measured 2D carrier densities is not com- 
pletely understood.) We conclude that existing studies 
are consistent with surface 2DEG and bulk superconduc- 
tivity in SrTiOs having a common origin. 

The strongly confined portion of the electron distribu- 
tion in the mid and high density regimes has significant 
{yz, zx} character. (See Figs. [7]&[9]). The increase in the 
density of states associated with these heavy 2D ba nds 
could account for the observed ferromagnetismJi^LUJ if it 
is describable by a Stoner criterion. The spatial separa- 
tion between the strongly confined and the low density 
tail portions of the 2DEG distribution may account for 
the coexistence of superconductivity and magnetism seen 
in some systems.^ This scenario should be compared 
with one in which superconductivity and magnetism both 



occur in strongly confined subbands; the two pictures 
should be experimentally distinguishable because of the 
strong influence of magnetism and spin-orbit coupling on 
superconducting properties^ in the spatially coincident 
case. 

The low density tail is a consequence of the property 
that the electric field vanishes far from the surface of 
a grounded system with a surface-bound 2DEG. The 
non-linear screening properties that we have discussed 
imply that a back-gate applied to the surface 2DEG 
to increase the strength of the electric field strength 
deep below the surface will have an exaggerated impact 
on the low-density tail of the distribution function and 
on the corresponding closely-spaced 2D subbands near 
the Fermi surface. A gate electric field with strength 
~ 10 14 cm~ 2 e oo /47re should essentially eliminate the tail 
region. Our prediction that superconductivity is associ- 
ated with the low density tail can therefore be tested by 
back gating which should suppress and eventually elimi- 
nate superconductivity^ without having a large influence 
on magnetism. Irrespective of the reliability of these pre- 
dictions, it seems clear that studies of the electronic prop- 
erties in dual-gated samples could be quite informative 
in building up a confident understanding of 2DEG prop- 
erties. 
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